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ABSTRACT 



The existence of black holes remains open to doubt until other conceivable 
options are excluded. With this motivation, we consider a model of a compact 
star in which most of the mass consists of dark particles of some kind, and a small 
fraction of the mass is in the form of ordinary nucleonic gas. The gas does not 
interact with the dark matter other than via gravity, but collects at the center 
as a separate fermionic fluid component. Depending on whether the dark mass 
is made of fermions or bosons, the objects may be called fermion-fermion stars 
or boson-fermion stars, respectively. For appropriate choices of the mass of the 
dark matter particles, these objects are viable models of black hole candidates 
in X-ray binaries. We consider models with a dark mass of 1OM and a range of 
gas mass from 10 _6 M Q to nearly 1M , and analyse the bursting properties of the 
models when they accrete gas. We show that all the models would experience 
thermonuclear Type I X-ray bursts at appropriate mass accretion rates. Since 
no Type I bursts have been reported from black hole candidates, the models are 
ruled out. The case for identifying black hole candidates in X-ray binaries as 
true black holes is thus strengthened. 



Subject headings: accretion — black hole physics — X-rays: binaries, bursts 
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1. Introduction 

Eighteen excellent black hole candidates have been discovered so far in X-ray binaries 
(XRBs, McClintock & Remillard 2003). The compact stars in these binaries have masses that 
exceed the maximum mass of a neutron star (Glendenning 2000); therefore, it is generally 
assumed that they must be black holes. The mass measurements certainly provide a strong 
argument for considering the objects to be black holes. However, it would be prudent to 
keep an open mind on this matter and to consider the possibility that the objects may be 
some kind of exotic stars that are composed of an as- yet unidentified form of exotic material. 
Until such a model is ruled out, the objects must be treated only as black hole candidates 
(BHCs), not true black holes. 

Glendenning (2000) has discussed a number of forms of exotic matter that might be 
present in compact stars, and describes models made up of these kinds of matter. All of 
the models have a maximum mass that is well under 3M . The models are, therefore, not 
relevant for BHCs in XRBs. There are other forms of compact stars, however, e.g. Q-stars 
(Miller, Shahbaz & Nolan 1998, and references therein), that are in principle consistent with 
the observed objects. Can one show that such objects are ruled out? One way to do this is 
to demonstrate that BHCs do not have hard surfaces. If we could show this, then the objects 
must have event horizons and must, therefore, be black holes. Beginning with the work of 
Narayan, Garcia & McClintock (1997, 2002) and Garcia et al. (2001), a number of studies 
have attempted to prove this result (e.g., Sunyaev & Revnivtsev 2000; Narayan & Heyl 2002; 
Done & Gierlinsky 2003). We are concerned in this paper with the work of Narayan & Heyl 
(2002), who argued that BHCs should exhibit Type I X-ray bursts if they have surfaces. 
We give an update on that work in §4 of this paper and show that the calculations rule out 
models of BHCs in which the objects have a hard surface on which accreting gas can collect. 

This still leaves open the possibility that BHCs may be made up of some kind of exotic 
dark matter with which normal gas does not interact. That is, the dark matter may be 
"porous" and permit accreting gas to fall through and to collect at the center. The dark 
matter and the fermionic gas would then behave as two independent fluids that interact 
only via gravity. Such models have been discussed in the literature, both for fermionic dark 
matter and bosonic dark matter (Lee & Pang 1987; Zhang 1988; Henriques et al. 1989, 
1990a,b; Jin & Zhang 1989, 1999). We refer to these objects as fermion-fermion stars and 
boson-fermion stars, respectively, where the first half of the name refers to the nature of 
the dark matter, and the second "fermion" in each name corresponds to the nucleonic gas 
component. If BHCs consist of either fermion-fermion or boson-fermion stars, would they 
have Type I bursts when they accrete gas from a normal binary companion star? We attempt 
to answer this question. 
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We describe in §§2,3 models of fermion-fermion and boson-fermion stars consisting of 
1OM of dark matter and various amounts of gas mass M 9 . We calculate how relevant 
properties of the gas sphere, such as the radius, surface gravity, redshift, etc., vary as a 
function of M 9 . Then, in §4, we describe the bursting properties of these models. We 
conclude in §5 with a discussion of the results. Appendix A relates the accretion efficiency 
of a fermion-fermion star to the redshift at the surface of the gas sphere. 



2. Fermion-Fermion Stars 

We consider a spherically-symmetric fermion-fermion star consisting of dark non-interacting 
fermions plus ordinary gas. We assume that the fermion particles have a mass rrif and that 
the star has a total of N* fermions, corresponding to a dark "baryonic mass" of M* = N*rrif. 
The ordinary gas component has a baryonic mass M 9 = N 9 m n , where m n = 936 MeV is 
the average mass of a nucleon, assuming the initial accreted gas contains 70% Hydrogen 
and 30% Helium and other elements (Narayan & Heyl 2003). The gravitational mass of the 
combined star is M grav . 

We solve for the equilibrium structure of the fermion-fermion star using the field equa- 
tions of General Relativity. In Schwarzschild coordinates, the line element of a spherically 
symmetric, time-independent spacetime can be written as 

ds 2 = -B{r)dt 2 + A{r)dr 2 + r 2 dtt 2 , (1) 

where A(r) = 1/(1 — 2GM(r)/r), G is the gravitational constant, and M(r) is the grav- 
itational mass enclosed within radius r. For an asymptotically flat spacetime, we have 
A(oo) — > 1 and B(oo) — > 1. For the specific case of the well-known Schwarzschild metric, 
B{r) = 1/A(r) = (1 — 2GM grav /r), where M grav is the total gravitational mass. Since our 
fermion-fermion star has a spatially distributed mass, we employ the more general metric 
written in equation (1). 

The metric functions A{r) and B(r) are determined by the Einstein field equations, 

R,u - \g, v R = 8nG[Tf u + T 9 U ], (2) 

where is the energy-momentum tensor of the fermions and T 9 V is that of the normal 
nuclear matter. We treat both the dark fermions and the gas as perfect fluids, 

T; V {P\ P l ) = {P l + P l )U,U v + g»„p\ (i = f,g), (3) 

where is the four- velocity of the fluid. From equations (1) and (3), we have the following 
structure equations, which are an obvious generalization of the Tolman-Oppenheimer-Volkoff 



-4- 



equations for the case of a single fluid, 

A'{r) = 8irGrA 2 (r)[p f (r) + p 9 {r)] - —\A(r) - 1], (4) 

r 

fl'(r) = 8nGrA(r)B(r)\p f (r)+p 9 (r)] + ^p-[A(r)-l], (5) 

j/(r) = -i|W[p/( r ) +p /( r )], (6) 

= -^|^[P 9 (0+P 9 W], (7) 

JV /# (r) = 47rr 2 v / A(r)n / (r), (8) 

N 9 '(r) = 4vrr 2 v / Z(r)n 9 (r), (9) 



where primes denote derivatives with respect to radius r, p-'(r) and p 9 (r) are the proper 
mass densities of the fermions and the gas, p* (r) and p 9 (r) are the corresponding pressures, 
and N*{r) and N 9 (r) are the total numbers of fermions and gas nucleons enclosed within 
radius r. 

The boundary conditions are as follows, 

A(0) = 1, 5(0) = 6o, /(0) = p{, p 9 (0) = p 9 c , N f (0) = 0, N 9 (0) = 0, (10) 

where p[ and p 9 are the central pressures of the fermion fluid and the gas, respectively. 
Because the equations are linear in B(r), b can be assigned an arbitrary value and later 
rescaled so as to satisfy the asymptotic flatness condition, B(oo) = 1. The two central 
pressures are free parameters that are determined by applying boundary conditions at large 
r on the total baryonic masses of the two fluids. 

In order to solve the above coupled equations, we need to know the equations of state 
of the fermions and the gas, 

P* = P>*), (11) 
P l =p\n l ), (12) 

where n % is the number density of particles of type % (= /, g). For simplicity, we treat the 
dark fermions as a non-interacting gas whose equation of state (EOS) may be written in 
parametric form (see Chandrasekhar 1935; Oppenheimer & Volkoff 1939), 

4 

rrir 

P = ^( sinht - t )' ( 13 ) 

l TTl f 

P = -^(sinht-8sinh(t/2) + 3t), (14) 
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where the parameter t is defined by 

1/2' 



/ = I. log < ^ + 

m f 



1 + 



„ 2 

Pf 
m f 



(15) 



and p_F is the Fermi momentum. The number density of fermions is 



n< = A.. (16) 

For the normal gas component, we use the FPS equation of state (EOS, Lorenz, Ravenhall 
& Pethick 1993) when the density is below the nuclear density p nuc , matching it smoothly to 
the relativistic mean field theory (RMFT) EOS at higher density, taking the so-called GM1 
set of the coupling constants (Glendenning & Moszkowski 1991). In practice, the maximum 
central density of the gas in the models is below 2p nuc , so different sets of the coupling 
constants give almost the same results. 

Starting from r = 0, the equations for the fermion pressure and fermion number are 
integrated out until the pressure goes to zero. The radius at which this happens is identified 
as the surface of the fermion star. So too for the gas. In calculating the fermion-fermion 
models, we first assumed that we have a pure dark fermion star with no gas. We adjusted the 
central pressure pf(0) of this model such that the asymptotic gravitational mass M^ rav (r — > 
oo) is equal to 10M Q . We selected the mass of the fermion particles to be 223 MeV. For this 
choice of mass, the maximum mass of a pure fermion star is 12.61M . A pure fermion star 
with a mass of 10M Q made of these particles is fairly compact: Rf = 252 km, Rf/Rs = 8.56 
in Schwarzschild units. Therefore, the model is interesting as a possible description of a 
black hole candidate in an X-ray binary. Figure 1 shows the metric quantities A(r) and B(r) 
as functions of radius for this model. 

We then added different amounts of gas to the above model, from baryonic mass M 9 = 
1O~ 6 M upto the maximum allowed mass beyond which the object becomes a black hole, and 
calculated the structure of the combined fermion-fermion star. In each model, we adjusted 
the central pressures p f and p 9 such that the total number of dark fermions is equal to the 
same obtained for the pure fermion model, and the baryonic mass of the gas is equal to 
the desired value. Figure 1 shows the results for the specific case when the baryonic mass 
of the gas is equal to 0.7 M Q . Note that the gas sphere is very much more compact than 
the fermion sphere in which it is embedded. This is the case for all the models we have 
calculated. 



For each choice of the baryonic mass of the gas M 9 , we obtained the radius R g of the 
gas component (the radius at which the gas pressure vanishes), the surface gravitational 
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acceleration g, and the redshift at the gas surface z. The redshift z(r) at radius r is defined 



as 



z(r) = -j== ~ I- (17) 
For g, we rewrite equations (6) and (26) in the following form, 

dp 9 l-B'(r) 1 / p 9 
~d£ " ~2B(r) ^/Mr) \ + P 9 

which defines the effective gravitational acceleration. At the surface of the accreted gas, 
particles are nonrelativistic, thus p 9 / p 9 <C I. The quantity £ in equation (18) is the surface 
mass density which is defined by dE(r) = p?{r) y/A(r)dr. We also calculated an effective 
accretion efficiency rj = L acc /Mc 2 , which measures the fraction of the rest mass energy of 
accreting matter that is released when gas accretes on the fermion-fermion star. This is 
obtained in terms of the gravitational mass of the star as follows: 

dM9 v ' 

We show in Appendix A that rj is given very simply in terms of the redshift at the surface of 
the gas sphere. The relation we derive is a relativistic generalization for a two-fluid system 
of a Newtonian result obtained by Rosenbluth et al. (1973). 

Figure 2 shows the results of the calculations. The four panels show the variations of 
the radius of the gas sphere, the gravitational acceleration at its surface, the redshift, and 
the accretion efficiency, as functions of the gas baryonic mass. With increasing gas mass, the 
radius of the gas sphere decreases, the gravitational acceleration increases, and the redshift 
and accretion efficiency increase. The wiggle in the curves at small masses corresponds to 
a switch between two solution branches. Note that the maximum mass in this sequence of 
models is O.782M ; beyond this mass, the object collapses to a black hole. 



3. Boson-Fermion Stars 

Boson stars are macroscopic quantum states resulting from the self-gravitational equi- 
librium of boson fields (see the reviews by Jetzer 1992; Lee 1992; Liddle & Madsen 1992; 
Mielke & Schunke 1998; Schunke & Mielke 2003). There is no concept of an equation of 
state for these systems, as they are pure quantum systems that are held up against gravita- 
tional collapse by the Heisenberg uncertainty principle. The pioneering studies in this field 
were done by Kaup (1968) and Ruffini & Bonazzola (1969). The main findings in these two 
seminal works are that (i) the mass of a boson star is of the order of Mpj/m b , and (ii) its 
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characteristic size is of the order of the de Broglie wavelength of bosons l/m b , where mj is the 
mass of the boson particle and Mp\ is the Planck mass. Colpi, Shapiro & Wasserman (1986) 
introduced the idea of self-interaction of the scalar particles and found that self-interacting 
boson stars have masses of the order of A 1 / 2 Mp l /mt ) , where A is a dimensionless quantity 
which characterizes the strength of the self interactions. For A 1 / 2 ^> 1, this scaling breaks 
down, and the mass is instead ~ Mj^/m 2 , which is similar to the result for fermion stars. 

Boson stars have several interesting/unique characteristics, including their transparency 
to photons and baryonic matter, their ability to be non-singular even for masses larger 
than the maximum mass of neutron stars (hence the interest in these objects as a model 
of BHCs), and the presence of a different metric than the Schwarzschild metric. Various 
authors have discussed observational consequences of boson star, such as the Cerenkov effect, 
gravitational lensing, the rotation curves of accreted matter, and the gravitational redshift of 
the radiation emitted within the effective radius of boson stars (Schunk & Mielke 2003). Also, 
the possibility of supermassive nonbaryonic stars (boson stars or neutrino balls) insteading 
of supermassive black holes existing in the nuclei of galaxies has attracted attention (Schunk 
& Liddle 1997, 1998; Torres, Capozziello & Lambiase 2000; Tsiklauri & Viollier 1998). It 
has been argued that the line profile of an emission line from an accretion disk around a 
supermassive boson star may have signatures that might help to identify the existence of 
boson stars (Lu & Torres 2003). 

The observational signatures of boson stars in X-ray binaries with BHCs has not been 
discussed very much in the literature. We show in this paper that such stars will produce 
Type I X-ray bursts with clear signatures. 

Henriques (1989, 1990a,b) considered the possibility of compact stars that contain both 
bosons and fermions and wrote down the structural equations for such objects. We borrow 
from their analysis in what follows. We analyse a boson-fermion star that consists of dark 
bosons plus ordinary fermionic gas. As in the case of the fermion-fermion star, we write the 
number of bosons as N b and the baryonic mass of the bosons as M b = N b mb- 

The energy-momentum tensor for bosons is completely different from that of fermions. 
For a massive self-interacting scalar field, the Lagrangian reads (Colpi et al. 1986) 



where A is the self-interaction coupling constant. According to Noether's theorem, the 
energy-momentum tensor can be written as 



L = ~ S TD lt 4>D v 4>-\ 



mU 2 ~ ~ A M\ 



(20) 



= D^D^ - -g^D^D^ + m 2 b 4> 2 + -A0 4 ). 



(21) 
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In order to solve the Einstein equations, it is reasonable to assume that 

<f)(r,t) = <&{r)e- iujt , 



(22) 



where <3>(r) is a real function. The structure equations can then be obtained from the Einstein 
equations (Henriques 1989), 



A'(x) = xA 2 (x) 



2p 9 (x) + 



Q 2 



B(x) 



l\a 2 + 4a 4 + 



a 



A(x) 



A(x) 



x 



[A{x)-1], (23) 



B\x) = xA(x)B(x) 



2p 9 {x) 



n 2 



B(x) 



A 



- 1 a z - -a* + 



a 



A(x) 



+ *M[A(s)-l](24) 



x 



p 9 (x) 
N b '(x) 
N 9 '(x) 



2 B\x) A'(x) 

° :;: "" { x + m^'2A(x) ]n 

lB'(x) 
~2 B(x] 

= \zx 2 t=. - ^A(x)a 2 , 



n 2 



B(x) 



- 1 \a - Aa 3 



.[p 9 (x)+p 9 (x)}, 



G y/Bl 



r 



-- A-kx 2 ^/ A(x)n 9 {x) ) 
where all physical quantities have been redefined in dimensionless form, 



(25) 

(26) 
(27) 
(28) 



x = rribr, a = \/4ttG$, Q = u/m b , A = A/ (AirGm 2 ,), 
p = AnGp/ml, p = AirGp/ml, n 9 = n 9 m^ 3 . 

The boundary conditions are as follows, 



(29) 



A(Q) 
p 9 (0) 



P 9 c 



B(0) = b , (7(0) = <7o, a'(0) = 0, 
N b (0) = 0, N 9 (0) = 0, n 2 (0)/b 



co- 



rn 



In addition, the total gravitational mass M grav is given by the asymptotic value of A at large 
r: A(oo) — > (1 — GM grav /r)~ 1 . As in the case of the fermion-fermion star, the equations 
are linear in b . Hence, this parameter may be given any value and normalized after the 
calculations have been done so as to achieve B(oo) = 1. Because the differential equation in 
a is second-order, we require two boundary conditions, one on a and another on a . As in 
the case of the fermion star, two parameters, namely o"o and p 9 , are determined by requiring 
the total boson and gas masses to equal specified values. In addition, there is an extra 
parameter c which behaves like an eigenvalue. This parameter is adjusted so as to obtain a 
wave-function a that is both nodeless and behaves regularly at infinity, i.e., a should remain 
positive for all x and should asymptote smoothly to as x — > oo. (If c has the wrong value, 
a either crosses and goes negative at a finite radius, or diverges as x — > oo.) 
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As in the case of the hybrid fermion star, we first solve for a pure boson star that has 
a gravitational mass of 1OM . We selected a boson particle mass of m b = 2.4 x 1CT 17 MeV 
and a coupling constant of A = 100. For these parameters, the maximum mass of a pure 
boson star is 12.57M Q . Figure 1 shows the variations of A(R) and B(R) with R for the pure 
boson star. The object does not have a precise surface since the wave-function extends to 
infinity. If we define the radius as the value of r at which 99.9% of the mass is enclosed, 
then the radius is 153.3 km. 

Having solved for the structure of the pure boson star, we then added various amounts 
of gas from 1O~ 6 M up to the maximum allowed mass of 0.863M Q (above which the object 
becomes a black hole) and solved for the structure of the combined boson-fermion star. 
Figure 1 shows the solution for the metric functions when M 9 = 0.7 M Q . Note that, even 
though the boson component does not have a precise surface, the gas component behaves 
like a fermion object and does have a surface at which the pressure goes to zero. The radius 
of this surface, the surface gravitational acceleration, the redshift and the accretion efficiency 
are shown as functions of M 9 in Fig. 3. 



4. Type I X-ray Bursts in Black Hole Candidates 

If black hole candidates (BHCs) in X-ray binaries are not true black holes with event 
horizons, then they must have surfaces of some kind. As discussed in §1, we may consider 
two general classes of models. In one class of models, BHCs are made of exotic matter and 
have a hard surface at some radius -Rbho When gas accretes on such a star, it is arrested 
at the surface and piles up as a layer of normal matter. Narayan & Heyl (2002) showed that 
the accreting layer will generally produce Type I X-ray bursts. 

Type I bursts in accreting neutron stars were originally discovered by Grindlay et al. 
(1976) and have been studied intensively for many years (see Lewin, van Paradijs & Taam 
1993 for a review of the observations, and Bildsten 1998 for a discussion of the physics). 
When gas accretes and accumulates on a compact star and is compressed by the strong 
gravity, thermonuclear reactions are ignited in the gas. In many cases (but not all, see 
Narayan 2003) the nuclear reactions are unstable and all the accumulated nuclear fuel is 
burned explosively within a short time. The result is a flash of thermal X-ray emission from 
the surface of the star, which is called a Type I X-ray burst. The theory of these bursts has 
been developed by a number of authors, e.g., Hansen & van Horn (1975), Woosley & Taam 
(1976), Joss (1977), Fujimoto, Hanawa & Miyaji (1981), Paczyhski (1983), Fushiki & Lamb 
(1987), Taam et al. (1993), Cumming & Bildsten (2000), Zingale et al. (2001), Narayan & 
Heyl (2003). The work described here is based on the method developed by Narayan & Heyl 
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(2003), which has some advantages over the earlier work. 

Figure 4 shows the results for the case of a 1OM BHC with an interior temperature of 
10 7 - 5 K (a reasonable choice, see Narayan & Heyl 2002). We consider different choices for 
the radius -Rbhc of the black hole candidate, from -Rbhc = (9/8)i?s, the smallest allowed 
value for a model in which there is no density inversion with increasing radius (Weinberg 
1972), up to -Rbhc ~ 3-R#, where Rs is the Schwarzschild radius corresponding to the mass 
of the BHC. We have considered a range of mass accretion rates and determined whether or 
not the system will have bursts. The solid dots in Fig. 4 indicate regions of the parameter 
space where the system will produce bursts, and the empty regions indicate regions where 
the accreting gas burns steadily without producing bursts. It is clear that in this model 
bursts are very common over a wide range of accretion rates. 

Figure 5 shows the expected recurrence times of bursts and the likely durations of bursts 
(assuming the maximum luminosity during a burst is equal to the Eddington luminosity) for 
a typical model with -Rbhc = 2-Rs. It is clear that for a range of accretion rates from a tenth 
of Eddington upto almost the Eddington rate, bursts occur reasonably frequently and have 
substantial fluences. These bursts would be hard to miss. The fact that no bursts have been 
seen from any BHC system in this accretion luminosity range is thus significant and argues 
against such a model for BHCs (Narayan & Heyl 2002). For accretion rates below a tenth 
Eddington, the systems are still unstable (see Fig. 4), but the bursts are relatively rare (and 
have correspondingly very large fluences). The rarity of the bursts means that this region of 
parameter space may be less useful for setting constraints on models. 

Another possibility for BHCs is the case we have focused on in this paper, namely that 
the objects are fermion-fermion or boson-fermion stars in which the gas collects as a compact 
sphere at the center of a dark fermion or boson sphere. Will such objects have Type I bursts 
when they accrete? To answer this question, we have redone the burst calculations using 
the models described in §§2,3. The quantities required for the burst calculations are the 
ones shown in Figs. 2, 3, namely the radius of the gas sphere R g , the surface gravity g, the 
redshift z(R g ), and the accretion efficiency r). With this information, it is straightforward 
to take any given model of a fermion-fermion or boson-fermion star, assume a local surface 
mass accretion rate S (baryonic mass added per unit area per unit time as measured in the 
local frame at the surface of the gas sphere), and compute the thermonuclear stability of the 
accretion layer. 

The two top panels in Fig. 6 show the results for the fermion-fermion star models 
described in §2 and the two bottom panels correspond to the boson-fermion models described 
in §3. In each panel, the horizontal axis gives the accretion luminosity in Eddington units, 



where LEdd is defined at infinity for the total mass of the star: 




(31) 



K 



with k = 0.4 cm 2 g _1 (corresponding to electron scattering in fully ionized hydrogen). The 
vertical axis extends over the full range of gas mass that we have considered in the hybrid 
models. The symbols (filled and open circles) in the panels correspond to models that are 
predicted to exhibit bursts and the empty spaces to models that are predicted not to have 
bursts. 

We see that there is a considerable range of accretion rate for which the systems will 
produce bursts. Therefore, the absence of bursts in BHCs argues against such a hybrid model 
for BHCs. A comparison of Figs. 4 and 6 shows that bursts occur in fermion-fermion and 
boson-fermion models for lower accretion rates compared to BHC models with hard surfaces. 
This difference has an interesting explanation. Recall that the Eddington luminosity was 
defined with respect to the total mass of the hybrid star, which is the obvious choice for an 
observer at infinity. However, the accretion and burst physics occur on the surface of the gas 
sphere which lies well within the dark matter sphere. Here, the local Eddington luminosity 
is given by LEdd, g as = ^R 2 g cg/ 'k, which is different from the Z^dd defined earlier. Redshifting 
this luminosity to the observer at infinity gives 



which is the physically correct luminosity to use as the scale when analysing bursts from 
these objects. We find that all the models we have analysed are unstable for all accretion 
rates up to the local Eddington limit. Such models are shown by the solid circles in Fig. 6. 
In fact, the models are generally unstable even for slightly super-Eddington local accretion 
rates, as shown by the open circles. 

Figure 7 shows the predicted burst recurrence times and burst durations for fermion- 
fermion and boson-fermion stars with a gas mass of O.3M (as an example). The recurrence 
times are a little long, but well within observational limits. The burst durations have been 
calculated with two different assumptions regarding the maximum luminosity during the 
burst. The solid circles assume that the maximum luminosity is equal to the full Eddington 
luminosity Z/Edd, while the open circles assume that the maximum luminosity is equal only to 
LEdd, g as,oo- Which is the appropriate one to use depends on whether or not the burst causes 
a substantial "radius expansion" of the radiating gas. This cannot be determined from the 
present calculations, but requires full time-dependent computations of the development of 
the burst. 



4nR 2 g cg 



(32) 



,gas,oo 



(1 + Z) 2 K 
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Regardless of the details, it appears that bursts should occur reasonably frequently on 
fermion-fermion and boson-fermion stars and that they should be easy to observe. The 
absence of these bursts means that such models may be ruled out for BHCs. 

5. Summary and Discussion 

For ordinary matter that is made up of baryons, and for any reasonable equation of 
state, the maximum mass of a compact degenerate star is about 3M©. Black hole candidates 
(BHCs) in X-ray binaries (XRBs) are more massive than this limit, so they cannot be 
baryonic. However, we cannot rule out the possibility that they are made up of some exotic 
non-baryonic form of matter. Until we can eliminate this possibility, we cannot be certain 
that BHCs are true black holes. We have considered two kinds of exotic stars in this paper. 

In one kind of model, we assume that the exotic star has a hard surface on which 
normal accreting gas can pile up. Such an object should behave very much like a neutron 
star except that it will have a different radius and surface gravity. Narayan & Heyl (2002) 
showed that such objects would exhibit Type I X-ray bursts when they accrete gas from a 
companion star. We have here repeated the calculations using the improved version of the 
burst calculations described in Narayan & Heyl (2003). The main freedom in the models is 
the radius that we assume for the surface of the BHC. We have tried a wide range (see Fig. 
4) and find that the results are not sensitive to the choice of radius. The calculations indicate 
that 10M Q BHCs with hard surfaces are unstable to bursts over an even wider range of mass 
accretion rates than neutron stars (compare Fig. 4 in this paper to Fig. 10 in Narayan & 
Heyl 2003). The burst durations and recurrence times for the BHCs are longer by a factor 
of a few than the corresponding time scales for neutron stars, but they are well within the 
reach of observations. 

An alternative variety of exotic star involves dark fermion or boson matter that does 
not interact with gas except via gravity. The lack of interaction means that accreting gas 
sinks to the center and forms a separate fluid component at the center of the dark matter 
distribution. The main free parameter in these models is the mass of the fermion or boson 
particles in the dark matter. There is not too much freedom, however, because we require 
the dark matter sphere to have a radius of no more than several Schwarzschild radii. This 
limit comes from the fact that observations of the X-ray emission of BH XRBs constrain the 
radius of the accretion disk to be quite small, no larger than about 6GM/c 2 (McClintock 
& Remillard 2003). High frequency quasi-periodic oscillations, in particular, require that 
most of the mass of the compact star should be inside this radius, which strongly constrains 
possible models of BHCs. In view of this constraint, we find that viable models are possible 
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only with fermions of mass ~ 200 MeV or bosons of mass ~ few x 10 -17 MeV. While we 
have done the calculations for one particular choice of the mass in each case (223 MeV and 
2.4 x 10~ 17 MeV, respectively), we expect the results to be similar for other choices within 
the relevant range. 

We have considered model objects consisting of 10M© of dark matter (fermions or 
bosons) and a wide range of mass in the gas component: 1O~ 6 M up to about 1M , the 
maximum mass beyond which the models become black holes. Over this entire range of 
models, we have calculated the thermonuclear stability of gas accreting on the surface of the 
gas sphere. In all cases, we find that Type I bursts are expected. Indeed, these models burst 
for all acretion rates up to the local Eddington rate as calculated at the surface of the gas 
sphere. (This is less than the Eddington rate for the combined gas plus dark matter object, 
see §4.) Moreover, the burst durations and recurrence times are not very unusual (Fig. 7 
shows results for the specific case when the gas mass is O.3M ), and the bursts should be 
easy to observe. 

What is a reasonable mass for the gas sphere in the models? The answer of course 
depends on how the particular exotic star is formed. If the star results from the death of 
a normal star, in which much of the gas is converted somehow to exotic matter, then one 
imagines that a reasonable fraction of the gas (perhaps several tenths of a solar mass?) 
might survive as normal matter. In this case, the upper range of gas masses that we have 
considered would be appropriate. On the other hand, if the exotic star is somehow born 
as a pristine purely dark matter object, then the only gas it would contain is whatever it 
accumulates through accretion. The typical mass transfer rate from the companion star in 
black hole XRBs is Mtransfcr ~ 10~ 10 - 10~ 9 A^yr" 1 (King, Kolb & Burderi 1996), and the 
typical lifetime of the systems is ~ 10 8 — 10 9 yr. Thus, we expect the accreted mass to be 
~ O.IMq. This is again near the upper end of the range we have considered (we went as low 
as 1O~ 6 M , see Fig. 6). 

The basic result of our study is that, for any kind of exotic star that we can consider, 
Type I X-ray bursts are expected to be present. Can these bursts be seen, and if so why 
have they not been observed? Tomsick et al. (2003) recently reported detecting bright 
~ 100 s flares in the X-ray flux of the BH XRB XTE J1650-500. These particular flares 
had nonthermal energy spectra, so they were not Type I X-ray bursts (which always have 
thermal spectra as one would expect for emission from a stellar surface). However, this work 
demonstrates that bursts are easy to detect in BHCs. The fact that no bona fide Type I 
burst has been reported in any of the dozen or more BHCs is thus highly significant. It 
would be useful to go back to the archival data on these systems and to derive quantitative 
limits on the burst rate. 



-14- 



To our knowledge, only one study so far — Tournear et al. (2003) — has specifically 
looked for bursts in BHCs. This work was done using data collected with the Unconventional 
Stellar Aspect Experiment and the Rosse X-Ray Timing Explorer. The authors derive a 95% 
upper limit on the burst rate in BHCs of 2.0 x 10 -6 bursts s -1 , compared to a mean detection 
rate in neutron star systems of (4.1 ±0.8) x 10 -5 bursts s -1 , which is about 20 times higher. 
Since our calculations indicate that the burst rate in BHCs (if the objects have surfaces) 
should be only a factor of a few lower than in neutron stars, the upper limit that Tournear 
et al. have obtained is already quite interesting. A longer observation might lead to quite 
strict limits on bursts in BHCs. 

The most obvious explanation for the lack of bursts in BHCs is that the objects are true 
black holes, with event horizons. Narayan (2003) has discussed various other explanations for 
the lack of bursts in BHCs and argues against all except the following two. One possibility is 
that BHCs are exotic stars with non-interacting dark matter through which gas can percolate 
down to the center. This is the model that we have analysed in the present paper. We have 
shown that such a model will produce observable bursts and is therefore not viable. The 
other possibility discussed by Narayan (2003) is that the unstable accretion layer on the 
surface of a BHC does not burst coherently in large explosions, in which the entire surface 
participates, but rather goes off in a series of random localized mini-explosions. The latter 
would be hard to distinguish observationally. Mini-explosions are expected if the burning 
front is unable to propagate rapidly over the surface of the star. It is hard to eliminate this 
possibility since the physics of deflagration fronts is not fully understood (but see Spitkowsky, 
Levin & Ushomirsky 2002). However, for a given nuclear fuel, we expect the propagation 
speed of the front to depend primarily on the local surface gravitational acceleration. Since 
g on the surface of a BHC is somewhat below that on a neutron star, and since both neutron 
stars and white dwarfs (which have much smaller values of g) are able to burst coherently, 
it is hard to see why BHCs should have any particular difficulty in producing large coherent 
bursts. 

Another caveat worth mentioning is that, in Fig. 4, we did not consider radii for the 
BHC smaller than (9/8)Rs- This radius is the smallest allowed within General Relativity for 
an object with a physical density distribution that satisfies dp/dr < (Weinberg 1972). If we 
allow a positive density gradient — which risks the Rayleigh- Taylor instability and requires 
a non-monotonic behavior of pressure with density — then smaller radii are possible. For a 
sufficiently small radius, the redshift would become very large and it would be possible to 
hide bursts from the view of the observer. As Abramowicz, Kluzniak & Lasota (2002) have 
noted, the gravitational condensate star described by Mazur & Mottola (2002) is one such 
model that can have an extremely large redshift. However, this particular model has serious 
conceptual problems. 
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Finally, we have assumed that the gas is fully degenerate, whereas one might wonder 
whether the gas component could remain hot and extended, and be thermally supported, for 
the entire age of the system (like low-mass brown dwarfs which retain their internal energy 
for nearly a Hubble time). Fortunately, this kind of a model may be ruled out. As mentioned 
earlier, X-ray observations indicate a small radius for the accretion disk, which means that 
the gas sphere must have an even smaller radius. For such radii, the virial temperature is 
very high (~ 10 11 K), and the density of the gas is also large. Neutrino cooling is then highly 
effective and the gas loses most of its thermal energy in a very short time (see Shapiro & 
Teukolsky 1983 for a discussion of the related problem of neutron star cooling). The gas 
thus quickly achieves degenerate conditions, where our calculations apply. 

The authors thank Jeremy Heyl for useful discussions. This work was supported in part 
by NASA grant NAG 5-10780 and NSF grant AST 0307433. Y.F.Y is partially supported 
by the Special Funds for Major State Research Projects, and the National Natural Science 
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A. Accretion Efficiency of a Fermion-Fermion Star 

We begin by quoting a result from Rosenbluth et al. (1973). Consider a gravitating star 
that is made of a single fluid whose internal energy e is a function only of the density p and 
satisfies de = —pd(l/p). If E(M) is the binding energy for a given mass M, Rosenbluth et 
al. (1973) show via a Newtonian analysis that 

= ®s, (Al) 

where $ s is the gravitational potential at the surface of the star. Thus, if a small quantity of 
mass AM is added to the star, the change in energy is equal to $ S AM. Since this quantity is 
equal to the energy released in dropping the mass from infinity down to the surface, equation 
(Al) implies that there is no additional internal energy released when the star adjusts to 
its new equilibrium. In a sense, the result is not surprising. Stars of different M all have 
the same specific entropy (because of the particular form of the equation of state), so there 
is neither an excess nor a deficit of specific entropy in going from one stellar mass M to 
another. 

We are interested in generalizing the above result for the relativistic problem, and also 
for the case of two fluids. Let us first consider a general relativistic star with a single fluid. 
Let M grav and Mb ar be the gravitational and baryonic masses of the star and let R be its 
radius. Let the fluid consist of particles of mass m with a number density n(r). Write the 
metric function A{r) as 



A(r) = 



2M(r) 



A(R)=[l- 2 -^) , (A2) 



This preprint was prepared with the AAS IATgX macros v5.2. 



-19- 



where we have used gravitational units with G = c = 1. The gravitational and baryonic 
masses of the star are given by 

f R 

M grav = / 4irr 2 pdr, (A3) 
Jo 

r R 

M bar = / Aixr 2 mn(l-2M/r)- l/2 dr. (A4) 
Jo 

In equilibrium, for a given baryonic mass, the star will take up the configuration with the 
lowest energy, i.e., the smallest value of M grav . Thus, for small variations of the equilibrium, 
we must have 

5M grav - A5M bar = 0, (A5) 

where A = <9M grav /<9Mb ar is a Lagrange multiplier. 

Equation (A5) must be valid for all first-order variations. Consider a specific variation 
consisting of a Dirac delta-function at the surface of the star, 

6n(r) = e6(r-R), Sp(r) = me5(r - R), (A6) 

where the second relation follows from the fact that n and p both go to zero at the surface, 
so that the fluid has no internal energy in this limit. Substituting (A6) into (A3) and (A4), 
we find 

5M grav = emAnR 2 , 5M hai = em4vr J R 2 (l - 2M grav / R)~ 1/2 . (A7) 
Substituting this in equation (A5), we then find 

where z(R) is the gravitational redshift at the surface of the star. This result is an obvious 
and natural generalization of the Newtonian result (Al). 

Consider now the problem of interest, a two fluid fermion-fermion star. Let us use 
subscripts 1 and 2 for the two fluids and assume that the outer radius R\ of fluid 1 is smaller 
than the radius i? 2 of fluid 2. (In the context of the main paper, fluid 1 is the gas and fluid 
2 is the dark matter.) Define the gravitational mass at a general radius r as 

M(r)= [ 4irr' 2 ( Pl + p 2 )dr', (A9) 
Jo 

so that the total gravitational mass of the star is 

M grav = M(R 2 ). (AlO) 
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The baryonic masses of the two fluids are given by 



Mi 



bar,l 



Mi 



bar,2 



/••Kl 

= / 47rr 2 m 1 n 1 (l-2M/r)" 1/2 rfr, 
Jo 

[Hi 

= / 47rr 2 m 2 n 2 (l-2M/r)- 1/2 rfr. 
Jo 



(AH) 



(A12) 



Requiring the total gravitational mass M grav to be a minimum for fixed baryonic masses of 
the two fluids gives the variational condition, 



<5M grav - A^Mb^ - A 2 <5M barj2 = 0, 

where Ai and A 2 are two Lagrange multipliers. 
Consider first a perturbation of the form 

Sni(r) = 0, Sn 2 (r) = e5(r — R 2 ). 



(A13) 



(A14) 



Substituting this in the above equations, and following the same steps as for the single fluid 
case, we directly obtain 



A 2 



I 2ilfg rav \ 



1/2 



R 2 J l + z{R 2 y 
which is identical to the result for a single fluid star. 
Consider next a perturbation 

6ni(r) = e5(r — Ri), 5n 2 (r) = 0. 



(A15) 



This gives 



5M(r) = em^TiRj, r > R u 
<5M grav = emiAirRl, 
5M harA = em x ^R\{\ - 2M/R l y 1 ' 2 , 

5M bar)2 = / 4vrrm 2 n 2 (l - 2M/ry 3/2 em 1 47rR 2 1 dr. 

Jri 



(A16) 

(A17) 
(A18) 
(A19) 

(A20) 



Substituting these relations in equation (A13) and making use of equation (A15), we obtain 
the following result for Ai: 



Ai 



2M(R 1 
Ri 



1/2 



1 - (l - ' 2 P ^rm 2 n 2 (l 



2M 



-3/2 



dr 



(A21) 
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This expression does not look as simple as equations (A8) or (A15), but in fact it simplifies 
considerably when written in terms of the redshift. 

To see this, start with the well-known Tolman-Oppenheimer-Volkoff structure equations 
for a relativisic star. Restricting ourselves to radii r > R u where there is only fluid 2, the 
equations take the form 



d[>2 
dr 



(p2 + P2)(M + 47rr 3 p 2 ) 
r 2 (l-2M/r) 

1 dp 2 



dr p 2 + p 2 dr 

Using the thermodynamic relation, 

P2 +P2 



dp2 

the above equations can be rewritten as 

n 2 



n 2 



dn 2l 



d 
dr 



P2 + P2 V r / 



-1/2- 



47rrn2 ( 1 — 



2M 



-3/2 



d\n\^B d ( ri2 
dr dr VP2+P2 
Integrating these equations and evaluating the result at r = Ri, we find 



n 2 (Ri)m 2 
P2(Ri) +P2OR1) 



2M(R 1 
Ri 



1/2 



2M a 



R2 



-1/2 



(A22) 
(A23) 

(A24) 



(A25) 
(A26) 



I I J 2Afg rav 



R2 



1/2 r R 2 



2M 



Vb(Ri) 



2M{R 2 ) 
R2 



1/2 



4:7irm 2 n 2 1 — 
Ri V 

n 2 (R 1 )m 2 



-3/2 



p 2 (i?i)+p 2 ( J Ri) 



, (A27) 
(A28) 



Finally, substituting (A27) in (A28) and comparing with (A21), we find that 



Ai 



grav 



<9M, 



bar,l 



Vb(Ri) 



i + z{R 1 y 



(A29) 



i.e., when written in terms of the redshift, the formula for Ai for the two-fluid case is identical 
to that for a single fluid. The accretion efficiency rj defined in the text is then given by 



7/ 



z(Ri) 



grav 



(A30) 
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We have confirmed equation (A29) by numerically computing the partial derivative 
9M grav /9M barjl and comparing it to 1/[1 + The agreement is very good. Indeed, 

numerical calculations show that equation (A29) is true also for a boson-fermion star, though 
we do not have a formal analytical proof for that case. 
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Fig. 1. — Variation of the metric potentials A{r) and B(r) with radius for representative 
models. Solid and dotted lines show the results for fermion stars and boson stars, respectively. 
The curves identified as "FF" show the results for fermion-fermion stars with 1OM of 
fermionic dark matter and 0.7 M & of normal gas, while the curves identified as "BF" show the 
results for boson-fermion stars with the same component masses. The curves not identified 
as FF or BF correspond to the pure dark matter stars with no gas. The squares show the 
surfaces of the gas and fermion components in the fermion stars, and the stars show the 
surfaces of the gas and boson components in the boson stars. For the boson component, the 
surface is defined as the radius that encloses 99.9% of the boson mass. 




Fig. 2. — Properties of the gas component in fermion-fermion stars. The four panels show 
the variations of the radius of the gas sphere R g , the surface gravitational acceleration g, the 
surface redshift z(R g ), and the dimensionless binding energy at the surface r], as functions of 
the baryonic gas mass M 9 . The fermionic dark matter component has a fixed mass of 10M Q 
in the models. 
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Fig. 3. — Properties of the gas component in boson- fermion stars. The four panels show the 
variations of the radius of the gas sphere R g , the surface gravitational acceleration g, the 
surface redshift z(R g ), and the dimensionless binding energy at the surface r], as functions of 
the baryonic gas mass M 9 . The bosonic dark matter component has a fixed mass of 10M Q 
in the models. 
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Fig. 4. — Burst stability of a 1OM black hole candidate with a hard surface. The accre- 
tion rate in Eddington units is shown along the abscissa and the radius of the object in 
Schwarzschild units is shown along the ordinate. The dots represent models that are un- 
stable to bursts, and the empty region (very close to the Eddington limit) corresponds to 
models that are stable and do not burst. Note that there is very little parameter space where 
bursts are absent. 
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Fig. 5. — Left: Shows the variation of the burst recurrence time t rec in days as a function of 
the accretion rate for a 10M Q object with a hard surface. A radius of 2Rs = 59 km has been 
assumed. Bursts are readily observable for accretion rates in the range from about O.lLEdd 
to about 0.7LEdd- Right: Shows the variation of the burst duration tburst in seconds for the 
same models. The durations have been computed assuming that the fluence in the burst 
emerges at the Eddington luminosity. 
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Fig. 6. — Top Left: Burst stability for fermion-fermion stars as a function of the accretion 
rate in Eddington units along the abscissa (calculated for the total mass of the object) and 
the mass of the gas component along the ordinate. The symbols correspond to models that 
are unstable to bursts, and the empty regions to models that are stable to bursts. Filled 
circles represent models that accrete at less than the local Eddington rate at the surface of 
the gas component. Open circles represent models that accrete at above the local Eddington 
rate but still well below the Eddington rate for the total mass. Top Right: Similar to Top 
Left, but for small gas masses (on a logarithmic scale) in the range 1(T 6 - 1O _1 M . Bottom 
Left, Bottom Right: Corresponding results for boson- fermion stars. 
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Fig. 7. — Top Left: Burst recurrence time t rec in days as a function of the accretion rate for 
a fermion-fermion star with a gas mass of M 9 = O.3M . Top Right: Burst duration ttmrst for 
the same models. The durations have been computed by assuming that the fluence in the 
burst is emitted at the Eddington luminosity; the filled circles are the results when the local 
Eddington luminosity at the surface of the gas sphere is used and the open circles when the 
Eddington luminosity of the total mass is used. Bottom Left, Bottom Right: Corresponding 
results for a boson-fermion star of the same mass. 



